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STEADY  STATE  TREATMENT  OF  RELATIVISTIC  ELECTRON  BEAM  EROSION 


I.  INTRODUCTION 

When  a  highly  relativistic  electron  beam  propagates  through 
initially  un- ionized  gas,  the  radial  forces  due  to  space  charge  and 
the  beam  magnetic  field  nearly  cancel  at  the  beam  head.  Consequently 
the  beam  head  is  not  self-pinched,  and  electrons  there  spread  out 
radially  at  a  rate  governed  by  the  beam  emittance.  Behind  the  freely 
expanding  head,  ionization  of  the  background  gas  raises  the  conduc¬ 
tivity  and  neutralizes  the  space  charge,  and  the  increasing  beam 
azimuthal  magnetic  field  pinches  the  beam  to  an  equilibrium  radius 
determined  by  the  balance  of  beam  pressure  and  the  radial  magnetic 
force.  When  the  background  gas  is  pre-ionized,  current  neutralization 
causes  a  similar  expansion  of  the  beam  head,  and  beam  pinching  then 
occurs  as  the  return  current  decays  resistively.  In  either  case  the 
beam  takes  on  a  characteristic  "trumpet"  shape. 

The  continual  expansion  of  the  beam  head  causes  the  ionization 
rate  near  the  head  to  drop.  Consequently,  the  region  where  the  pinch 
force  builds  up,  referred  to  here  as  the  "pinch  point",  moves  back 
toward  the  body  of  the  beam.  This  recession  of  the  pinch  point  from 
the  head  slows  down  as  the  head  radius  becomes  large,  until  a  second 
effect  becomes  important:  for  a  highly  relativistic  beam, energy  loss 
due  to  the  induced  axial  electric  field  ahead  of  the  pinch  point 
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reduces  the  relativistic  mass  of  beam  electrons  and  appreciably  increases 

radial  expansion.  After  an  initial  transient  period  of  head  expansion, 

the  two  effects  together  cause  an  almost  steady  movement  of  the  pinch 

point  back  into  the  beam  and  a  continual  effective  loss  of  the  head. 

This  beam  head  "erosion"  can  significantly  affect  the  stability  and 

propagation  characteristics  of  a  high-energy  electron  beam. 

Several  groups  have  studied  beam  head  erosion  and  have  informed  us 

1  2 

of  their  unpublished  results.  Both  Biegalski  and  Lee  have  derived 

analytic  estimates  of  the  erosion  rate  using  simple  energy-balance 

arguments.  To  study  beam  structure  and  verify  the  predicted  erosion 
2  3 

rate  scaling,  Lee  and  Brueckner  have  made  time-dependent  calculations 
of  axial  beam  profiles  using  radially-averaged  beam  equations.  Similar 

4 

calculations  have  been  carried  out  by  Hubbard  et  al.  and  by  Johnston 

5 

et  al.  with  more  elaborate  models  of  conductivity  production  and  the 
beam  fields. 

The  present  work  obtains  erosion  rates  for  relativistic  electron 
beams  by  a  different  method.  Since  the  erosion  rate  reaches  an  approx¬ 
imately  constant  value  when  the  head  has  expanded  sufficiently,  we 
assume  that  the  beam  axial  profile  is  stationary  in  a  frame  of  reference 
moving  with  the  pinch  point.  This  assumption  leads  to  a  time- 
independent  set  of  radially-averaged  beam  equations  that  we  solve 
numerically  to  obtain  both  the  pinch  frame  velocity  and  the  axial 
beam  profile.  This  method  requires  substantially  less  computation  than 
time-dependent  profile  calculations  and  gives  the  erosion  rate  un¬ 
ambiguously.  Details  of  the  model  are  presented  in  Sec.  IZ  along  with 


the  methods  used  to  solve  the  beam  equations.  Section  ZZI  gives 
results  of  the  model:  An  analytic  approximation  of  the  beam  head  erosion 
rate  is  derived  and  used  to  reduce  the  beam  equations  to  a  simpler 
dimensionless  form.  These  equations  show  that  the  beam  axial  profile 
depends  principally  on  the  dimensionless  parameters  that  determine  the 
gas  conductivity.  We  propose  analytic  scaling  laws  for  the  beam 
radius,  pinch  force,  and  conductivity,  and  compare  the  analytic  pre¬ 
dictions  with  extensive  numerical  solutions.  Section  ZV  then  summarizes 
the  conclusions. 
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II.  MODEL 


A.  Assumptions 

The  configuration  considered  here  is  an  axi symmetric  beam  of 

relativistic  electrons  moving  paraxial ly  in  the  positive  z  direction. 

For  convenience,  the  beam  head  position  at  time  t  -  0  is  chosen  to  be 

z^O,  and  the  beam  source  is  at  z  =  -°°.  There  are  no  externally- 

generated  fields.  We  take  the  beam  current  I.  to  be  small  compared 

b 

with  the  Alfven  limit  IR  =  17  (y2  -  l)1/2  kA,  where  y  is  the  relativ¬ 
istic  factor  of  typical  beam  electrons.  Since  this  requirement  ensures 
that  the  transverse  velocity  of  beam  particles  vA  is  small  compared 
with  the  axial  velocity  v^ ,  we  make  the  "paraxial  beam”  approximation , 

consistently  neglecting  terms  of  the  order  of  v2 /v2  compared  with  unity. 

i.  z 

Also,  we  assume  that  y  i  5  everywhere,  so  that  v^a  c,  and  the  relative 

axial  motion  of  beam  electrons  is  negligible. 

In  this  paper  we  do  not  consider  instabilities,  such  as  the 

well  known  resistive  hose  instability.  The  development  of  a  beam  is 

then  a  time-dependent  axisyrametric  process  involving  both  the  radial 

and  axial  dimensions.  In  this  work  we  make  three  approximations  that 

allow  a  one-dimensional  time- independent  formulation: 

1.  We  take  I.  ,  y,  and  the  beam  emittance  to  be  time- independent 
o 

at  the  source  throughout  the  pulse,  and  we  neglect  changes  in  beam 
emittance  and  energy  due  to  collisions  and  radiation.  These  single¬ 
particle  loss  processes  are  instead  assumed  to  cause  only  gradual  para¬ 
metric  changes  in  the  beam  erosion  rate  and  axial  profile.  However, 
changes  in  y  due  to  the  collective  electromagnetic  fields  of  the  beam 
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are  calculated  self-consistently;  these  effects  are  very  strong  near 

the  pinch  point  and  are  essential  features  of  bean  head  evolution. 

2.  All  bean  quantities  are  assumed  to  depend  on  z  and  t  only 

through  the  combination  £  =  V  t  -  z,  where  V  =  Be  <  v  is  the  effec- 

P  P  z 

tive  pinch-point  velocity  in  the  laboratory  frame  of  reference.  This 

assumption  implies  that  the  beam  structure  is  stationary  in  a  frame 

moving  with  velocity  V  .  Although  the  beam  head*  located  as  £  ,  is 

P  o 

obviously  moving  relative  to  this  pinch  frame  when  V  <  v  »  we  show 

p  z 

later  that  the  expanded  beam  head  near  £q  has  little  effect  on  the 

pinch  structure  when  the  radius  at  £  has  become  large  compared  with 

o 

the  pinched  radius.  Also,  we  assume  that  the  beam  tail,  at  £T =  £  +  L, 

L.  o 

is  well  behind  the  pinch  point.  Mathematically,  the  stationary- 
profile  assumption  implies  that 


9  Q  9 

at  *  Bc  9£ 

3 _ 3_ 

3z  3£  ' 

where  B  =  V  /c  must  still  be  determined.  Using  the  fact  that 
P 

v  %  c  »  v.  for  beam  electrons,  we  also  write  all  total  time 

z  * 

derivatives  as 


d_ 

dt 


9  *  3 

at  c  az  = 


(2) 


3.  The  beam  radial  profile  is  assumed  to  change  self- 
similarly  as  the  beam  radius  expands.  We  therefore  model  the  beam 
fields  and  current  density  and  the  background  conductivity  by 
functions  of  the  form 
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flr.Z )  -  f(0,O  F  (r/R) 


(3) 


where  r  is  the  radial  coordinate,  R(€)  is  the  mean  beam  radius ,  and  the 
radial  envelopes  F(r/R)  are  functions  with  F(0)  =  1.  This  assumption 
avoids  a  self-consistent  calculation  of  the  beam  radial  structure  and 
permits  analytic  radial  averaging  of  the  beam  equations.  We  find  that 
calculated  head  erosion  rates  and  axial  profiles  are  insensitive  to  the 
choice  of  radial  dependences. 

with  these  assumptions,  we  can  write  the  equations  for  the 
beam  radius  and  energy  and  the  radially-averaged  fields  and  conductivity 
as  a  set  of  coupled  ordinary  differential  equations  in  the  variable  Z , 
with  the  dimensionless  pinch-point  velocity  B  appearing  as  a  "nonlinear 
eigenvalue". 
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To  clarify  the  assumptions  involved,  we  review  the  derivation 
of  the  equations  for  the  beam  radius  and  mean  particle  energy  from  the 
single  particle  equation  of  motion.  For  a  collisionless  nonradiating 
electron  in  electric  field  E  and  magnetic  field  B,  the  relativistic 
equation  of  motion  is 


St  w 


(v  x  B  \ 

£  +  V1)' 


(4) 


where,  -e  and  m  are  respectively  the  electron  charge  and  mass,  and 
Y  =  (1  -  y y/c2)”1/2  is  the  ratio  of  total  particle  energy  to  rest 
energy  me2 .  The  dot  product  of  y  with  (4)  gives  the  energy  equation 


dy 

dt 


(5) 


If  terms  of  the  order  of  (v  /c)2  and  y-2  are  discarded,  then  (4)  yields 


d_ 

dt 


(6) 


where  to2  =  -  e (Bn -E  )/(mr)  is  the  nonrelativistic  betatron  frequency, 
p  or 

and  (5)  reduces  to 


a? _ e__ 

dt  ~  me  z  ’ 


(7) 


Since  beam  electrons  traverse  much  of  the  beam  cross-section 
while  performing  betatron  oscillations,  it  is  a  reasonable  approxi¬ 
mation  to  neglect  any  correlation  between  the  energy  mc2Y  of  an  electron 
and  its  radial  position  r.  (The  limitations  of  this  approximation,  and 
effects  due  to  spread  y  and  y-x  correlation,  are  discussed  in  Ref.  6.) 


We  therefore  assume  y  is  a  function  of?  only  and  write  an  averaged  form 


*L 

dt 


(8) 


where  &  is  the  average  axial  electric  field  acting  on  electrons  in  a 
thin  beam  segment.  Explicitly,  we  define  the  average  of  any  quantity 
f(r,C)  over  the  beam  cross  section  as 


<f>"  %  \  **  r  nfaf' 
0 


(9) 


where  n^r,?)  is  the  beam  number  density,  and 

=  2tt  J*  dr  r  n^  (10) 

0 

is  the  number  of  beam  electrons  per  unit  length. 

To  derive  an  equation  for  the  root-mean-square  beam  radius  R, 
appropriate  moments  of  C6)  sure  averaged  over  the  beam  cross  section, 
and  self-similar  expansion  is  assumed  (but  it  is  not  necessary  to 
assume  any  particular  self-similar  form) .  The  resulting  envelope 
equation,  derived  in  more  general  form  by  Lee  and  Cooper7,  is 

it  (;r3  it  (*  f )  *  ^  v2  ]  -  °-  1111 

The  beam  is  assumed  to  be  in  quasi-static  force  balance  at  its  tail, 

S  *  ?  ,  so  that  the  time  derivatives  of  y  and  R  vanish  there.  Slow 

Jj 

variation  of  y  and  R  at  the  beam  tail  will  be  treated  in  Sec.  HE  as  a 
parametric  variation  of  the  boundary  conditions.  The  first  integral  of 
(11)  is  then 
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where  U  h  <r(B.  -  E  )>,  and  the  subscript  L  denotes  quantities  at  C  • 
®  *  L 

We  use  the  steady-state  relation  (2)  to  rewrite  (8)  and  (12) 


as  time- independent  equations  for  the  beam  radius  and  energy: 


(13) 


&  =  _ 2 _ 

d?  (1-B)mc2  cz 


(14) 


■ 


C.  Field  Equations 

Since  the  beam  magnetic  field  only  appears  in  (13)  combined 

with  the  radial  electric  field  in  U  =  <r(Ba-E  )>,  we  reduce  Maxwell's 

s  r 

equations  to  a  pair  of  equations  for  u  =  r(B.-E  )  and  E  .  In  the 

Sr  z 

absence  of  external  magnetic  fields  and  azimuthal  currents,  the  com¬ 
ponents  of  B  and  E  in  a  cylindrical  geometry  are  B  =  (0,B_,0)  and 
i-  =  (E^/OfE^),  and  the  pertinent  components  of  Maxwell's  equations  in 
Gaussian  units  are 


r  (rV  =  T~  (JZ  +  aE>  +  -  gT5- 
r  or  9  c  z  z  c  3t 


(15) 


!!§ 

3z 


.  3E 

—  (J  +  CTE  )  +  -  r— 
c  r  r  c  3t 


(16) 


3E 

_ i 

3z 


3E 
_ 2 

3r 


i!!l 

c  3t 


(17) 


where  J  and  J  are  components  of  the  beam  current  density  and  a  is  the 
r  z 

scalar  background  conductivity.  When  (15)  -  (17)  are  written  in  terms 
of  £  using  (1) ,  the  equations  can  be  combined  to  give 


(1  +  8)  (B.-E  )  +  ^2-  (B0-  E  )  -  —  J  +  r-2-  +  —  B 

3£  0  r  corcrSr  c9 


(18) 


**  3  3Bfl 

3 r  *  3?  (W  -  (1-0)  3T  {19) 


r  3r  (rB0) 


$ 


(20) 
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The  radial  beam  current  is  eliminated  by  invoking  charge 
conservation.  In  terms  of  £ ,  the  continuity  equation  is  written 


l  a 

-e8car  +  ?  V  -sr=0' 


and  with  as  -ecn^,  this  yields 


where 


j  =  £=§.  ii 
r  2irr  3£ 


I(r)  s  2ir  1  dr'r'j'  , 


and  J"  denotes  J  (r ") .  Using  (22),  we  rewrite  (18)  and  (19)  as  a 
z  z 

coupled  pair  of  equations  for  u  and  E 


ii*a\  3“  .  4irc  ..  2(1-B)  dl  .  8Ez  .  4tto  _ 

(1+B)  3c  +  ~T  u  ~  af  +  r  IT  +  T  rV  (23) 


-  |?  |  t«'-  a-eirsy]. 


with  BQ  given  explicitly  by  the  integral  of  (20) : 


2i+^f 

c  c  J 


dr'r'c'E'  +  8 
z 


tc  J  ar'r'Ei- 


Here,  E  is  taken  to  vanish  at  r  *  R  .  This  cutoff  radius  for  E  is 
z  c  z 

mathematically  necessary  to  give  the  beam  a  finite  inductance  per  unit 
length,  and  it  corresponds  physically  to  the  radius  at  which  the  con¬ 
ductivity  becomes  too  small  to  cause  charge  neutralization.  Normally 


R  is  determined  by  the  extent  of  induced  breakdown  around  the  beam 
c 
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head  and  is  large  compared  with  the  fully  pinched  radius  R^.  since 
the  cutoff  radius  enters  the  field  equations  only  logarithmically,  we 


take  R  to  be  constant.  When  1-S<<1,  (23),  (24),  and  (25)  reduce 
c 

approximately  to  the  single  equation 


If +  “  ["♦  It  j  **-*’<]  *  i  *  gi  f  r  ]  •  (26) 


with  the  supplementary  relation 


(27) 


For  typical  high-y  beams,  we  find  that  1-8  <  0.2,  justifying  the  use  of 
(26)  and  (27).  We  note. that  1-8  is  retained  in  (13)  and  (14),  where 
it  is  obviously  important. 

To  obtain  equations  for  0  2  <u>  and  €  ,  we  choose  specific 

z 

self-similar  expressions  for  n^,  u,  and  c  with  the  form  of  (3),  and 
average  the  field  equations  (26)  and  (27)  according  to  (9)  ..  since  the 
number  density  of  a  relativistic  beam  in  an  initially  uniform  neutral 
gas  changes  from  an  approximately  Gaussian  radial  envelope 

Vr'°  "  ~r2To  exp  (-r2/R2)  (28) 

0 

in  the  freely  expanding  head  to  a  Bennett  profile 


Nb  1 _ 

nb(r,C)  -  7tr2'(s)  (l  +  rVR^Z 


in  the  pinched  region,  either  of  these  forms  might  be  used  for  n^.  In 
most  numerical  calculations,  we  use  the  Gaussian  profile  (28).  For 
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comparison,  calculations  are  also  carried  out  using  a  cutoff  Bennett 
distribution , 


1.255 


r  «  1.98R, 


n^r,?)  *  0, 


r  >  1.98R, 


(30) 


and  a  square  distribution, 

Nb 

Vr'5)  *  ' 

n^r,?)  *  0, 


r  «  21/2R, 
r  >  21'^2R, 


(31) 


where  the  cutoff  radii  have  been  chosen  so  that  <r2>  *  R2.  We  show 
later  that  the  calculated  erosion  rates  and  axial  profile  are  not 
sensitive  to  the  form  chosen  for  n.  » 

D 

The  pinch  force  u  is  assumed  to  have  the  same  radial  dependence 
as  rB  for  the  beam  current  alone.  This  choice  is  equivalent  to 
assuming  that  fractional  neutralization  of  the  beam  current  and  space 
charge  is  uniform  over  the  beam  cross  section.  For  a  Gaussian  beam 
density  profile,  u  is  then  given  by 

u  =  2U[1  -  exp  (-r^/R2) ] ,  (32) 


and  analogous  analytic  expressions  are  readily  found  for  the  Bennett 
and  square  beam  profiles. 

For  a  beam  propagating  into  an  initially  un-ionized  or  weakly 
ionized  gas,  the  gas  conductivity  is  due  primarily  to  direct  ionization 
by  the  beam.  The  radial  conductivity  profile  is  therefore  expected  to 
be  qualitatively  similar  to  the  beam  density  profile.  However,  the 
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conductivity  profile  at  any  c  is  typically  broader  than  the  beam 
density  profile  at  the  same  position,  because  the  gas  ionization  is  due 
to  the  more  expanded  part  of  the  beam  ahead  of  5.  In  addition,  such 
effects  as  electric  breakdown  and  beam-generated  radiation  can  result 
in  broad,  weakly-conducting  radial  wings.  As  with  the  beam  density 
profile,  we  do  not  attempt  a  self-consistent  calculation  of  the  radial 
conductivity  structure.  Instead,  we  specify  an  analytic  self-similar 
form  for  the  radial  dependence  of  a  and  examine  the  sensitivity  of 
calculated  results  to  this  choice.  For  the  Gaussian  density  profile 
(29) ,  we  assume  the  form 

o (r,0  *  a (0,5)  exp  t-ar2/*2 (?) ] ,  (33a) 

where  a  is  a  constant  between  zero  and  unity  which  can  be  used  to  vary 
the  breadth  of  the  conductivity  profile  relative  to  the  beam  profile. 

We  find  that  the  results  of  numerical  calculations  are  quite  insensitive 
to  the  choice  of  a,  which  justifies  the  use  of  simple  self-similar 
models.  For  beams  with  the  cut-off  Bennett  density  profile  (30),  we 
use  the  unbroadened  conductivity  profile 

=  11 «)V  ■ for  r* 1-98"-  l33b> 

and  for  the  square  density  profiles  (31)  we  take 

a(r,0  *  a(0,O,  for  r  21/2r  (33c) 

The  calculation  of  a (0,0  is  discussed  in  the  next  section. 

The  radially-averaged  quantity  k  =  4  ir<c>/c  occurs  frequently 
in  the  averaged  field  equations.  For  each  of  the  three  beam  profiles, 
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k(Z)  is  related  to  a  10,0  by 


<(0 

=  (l+ct)_14Trcr(0,O/c, 

Gaussian, 

(34a) 

ic(0 

=  (0.415)  4iro(0,O/c, 

cut-off  Bennett, 

(34b) 

<(0 

*  4ira  (0,O/c, 

square . 

(34c) 

With  these  assumed  radial  dependences,  a  straightforward  but 
lengthy  calculation  gives  the  radial  averages  of  (26)  and  (27),  For  each 
of  the  three  radial  density  profiles  the  equations  take  the  form 


d2U  dU  /*b  \ 

—  +  B  —  -  CTJ  =  -  U  j  , 


d  0  dR  ,  n  dU 

'  5?  <xo>  -  5  as  +  D  «  ' 


where  A  2  Zn(R  2/R2) .  For  a  Gaussian  density  profile  with  the 
c 

conductivity  profile  (34) ,  the  coefficients  in  (35)  are 

2  +  a  l\  .  »  A  +  a\  .  ..  "1 


rrsr  [x  *  *  V 


(37a) 

„ dR 
KRat  - 

(37b) 

n 

(37c) 

(37d) 

D  -  In  2  +  ye, 


where  y_  ■  0.5772  is  Euler's  constant.  For  the  cut-off  Bennett 
E 

density  profile  with  a  Bennett  conductivity  profile  of  the  same  width, 
we  find 
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(38b) 


B  -  1.59  +  (0.46X  +  0.15)k2R2  -  0.91kR  , 

aC 

C  *  <1.50  +  0.70k2R2)  |  ||  +  [°*45  J  fjT  _  °‘34  (i  If  )  ]  (38c) 

D  *  0.111.  (38d) 

When  the  density  and  conductivity  have  square  radial  profiles  extending 
to  21/2r,  the  coefficients  are 

A  ,  I  (X  +  |  -  ln2)KR2,  (39a) 

B-1+J-CX  +  I-  -  £n2)K2R2  -  I  kR  ?  ,  (39b) 

2  3  3  as 


(2+1  K2R21  I  §*  +  1  K  R  £* 

(2+3KR)  RdC+3K  V 


(39c) 


D  *  In  2  -  J 


(39d) 


.  Conductivity  Equations 


The  processes  governing  the  evolution  of  the  gas  conductivity 
can  be  modeled  at  various  levels  of  completeness  and  accuracy.  In  this 
work,  we  use  the  simplest  model  which  is  appropriate  for  weakly- ionized 
gases,  the  regime  of  principal  interest  to  us.  In  this  regime,  the 
conductivity  is  related  to  the  background  electron  number  density  ng  by 


a 


n  e*- 
e 

mv 

m 


(40) 


Here,  is  the  electron-molecule  momentum  transfer  collision  frequency, 
which  is  proportional  to  the  gas  density  and  only  weakly  dependent  on 
electron  temperature  T  .  For  dry  air, 

v  ss  1.4  x  1012  T  */2  p  sec  ^  ,  (41) 

m  e 

where  T  ,  in  units  of  eV,  is  between  0.5  and  10,  and  p  is  the  air 
e 

density  normalized  to  the  density  of  air  at  standard  temperature  and 
pressure * 

We  model  the  electron  density  as 


3n 

5T  '  V|jbl  '  VV  *  VTe’°V  (42> 

where  the  first  term  represents  direct  ionization  by  beam  electron 
impact,  the  second  term  represents  three-body  recombination  and  the 
third  term  represents  ionization  by  secondary  electron  intact,  the 
process  which  leads  to  electric  breakdown.  Since  is  independent  of 

A 

T  and  the  temperature  dependence  of  A?  is  fairly  weak,  we  take  both 
e  *■ 
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r 

! 


coefficients  to  be  constant.  In  contrast,  increases  sharply  with 

T  ,  particularly  for  T  2  5eV.  To  model  breakdown  accurately,  we 
e  e 

would  either  have  to  include  an  equation  for  Tg  in  the  formalism  or 

use  some  phenomenological  relation  giving  in  terms  of  quantities 

that  determine  Tg,  such  as  the  ratio  of  electric  field  to  gas  density. 

We  note,  however,  that  for  p  2  0,1,  secondary  ionization  is  usually 

only  a  small  correction  to  primary  beam  ionization,  except  in  the 

radial  wings  outside  the  beam.  The  reason  is  that  Tg  is  large  only  in 

front  of  the  pinch  point  where  substantial  E  and  E  fields  occur,  but 

ne  is  small  enough  here  that  primary  ionization  dominates  secondary 

production.  We  therefore  ignore  the  A3  term  and  use  temperature- 

independent  approximate  values  of  A.  and  v  .  We  then  use  (40)  and 

2  in 

(34)  to  rewrite  (42)  with  r  *  0  as  an  equation  for  the  radially 
averaged  quantity  k 


where 


xi*,bl 

R2 


2-C2  , 


4e2A^ 

\  9  nl  mc2v  /F  ' 
m 


m(v  /p)A 
m _ 

4ire2 


(43) 


(44) 


(45) 


and  and  n2  are  geometric  factors,  given  by 


n 


1 


(1  +  o)  Gaussian 
0.520  ,  cut-off  Bennett 


(46) 


0.5 


square 


II  +  cc  ,  Gaussian 

2.41  ,  cut-off  Bennett  (47) 

1  ,  square 

A  «■» 

Values  of  A,  ,  A, ,  and  v  /p  are  shown  in  Table  1  for  several  gases  at 
11  m 

T  =  2eV.  The  recombination  rate  has  little  effect  on  either  beam 
e 

head  erosion  or  beam  structure  near  the  head,  the  principal  subjects 
treated  here.  In  all  calculations  presented  in  Sec.  Ill  we  use  a 
typical  value  for  air  A^  =  8.7  x  10  We  note  that  as  £  -*■  <*>,  k 
approaches  an  asymptotic  value  given  by 


k  R 


00  00 


(48) 


Our  interest  is  principally  in  the  regime  where  tc^R^  »  1.  For  large 
?  the  beam  space  charge  is  then  neutralized,  the  beam  pinch  force  is 
magnetostatic,  and  the  magnetic  decay  length  is  long  compeared  to  the 
beam  radius.  For  air  at  standard  density,  <  R  >  100  if  1 1.  |  >  0.3kA, 

oo  oo  D 

and  the  assumption  that  the  air  be  weakly  (s  1%)  ionized  is  satisfied  if 
Ilj/R2  S  5  x  103  kA  cm'2. 
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.  Boundary  Conditions 


The  envelope  equation  (13) ,  energy  equation  (14) ,  field 

equations  (35)  -  (39) ,  and  conductivity  equation  (43)  form  a  complete 

set  of  equations  for  R,  y,  U,  £  ,  and  k.  Since  the  R  and  U  equations 

z 

are  second  order  and  the  £  equation  is  zero  order,  six  boundary 

z 

conditions  are  required.  For  a  highly  relativistic  beam,  the  fields 
ahead  of  ?  vanish,  so  appropriate  boundary  conditions  for  ic  and  U  are 


U(t  )  =  U  =  0 

(49) 

o  O 

<(C  )  =  <  >  0  , 

o  o 

(50) 

where  the  subscript  zero  hereafter  labels  quantities  at  £  *  £  .  The 

o 

initial  conductivity  may  be  chosen  to  model  any  degree  of  pre¬ 
ionization. 

To  obtain  boundary  conditions  at  the  beam  tail  t  ,  we  assume 

Zj 

that  the  tail  is  far  enough  behind  the  pinch  point  that  the  beam  is  in 
quasi-static  force  balance  at  £  .  The  appropriate  boundary  conditions 
are  then 


y«l)  =  v 


R(?l)  =  V 


dR 

dC 


_  r2  §_  1 

d5  (yU)1/z 


U(CL) 


u. 


w  ' 


(51) 

(52) 

(53) 

(54) 


where  e 


Vl*L2 


is  the  beam  emittance  at  £  . 

L 


To  relate  e  to  the 
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beam  source  characteristics,  we  assume  that  electrons  at  C  have  not 

L 

undergone  significant  ohmic  energy  loss.  This  assumption  is  usually 
a  good  approximation  because  ohmic  dissipation  is  found  to  be  concen¬ 
trated  ahead  of  the  pinch  point.  If  single  particle  energy  loss 

mechanisms  and  beam-gas  scattering  are  then  neglected,  y  and  t  are 

Xj 

equated  with  the  values  at  the  beam  source.  If  energy  loss 

by  ionizing  collisions  and  bremsstrahlung  are  included,  the  slow 

variation  of  y  with  the  distance  of  the  beam  tail  from  the  source  is 
9 

given  by 


dy  y  /  e2  \2  -  /nmc2y3/'2\ 

&  ■  -  >r  -4"’02pr )  u ^  w—  ) 


(55) 


and  the  slow  variation  of  e,  due  to  scattering,  is  given  by 


7,10 


de2  4itR2c2  he 


dz 


where  ^  is  the  electron  bremsstrahlung  radiation  length,  given  by 


(56) 

9 


X  _1  =  4  n  Z 
r  o 


2/2  \  2 

(Z  +  1)  r-(^T  I  (183Z-1/3),  (57) 

nc  ymc4  J 


Wr«13eV,  and  nQ  and  Z  are  respectively  the  number  density  and  effective 

nuclear  charge  of  the  background  gas.  In  air  at  standard  density  X^«  270 

meters,  and  the  second  term  in  (55)  is  of  the  order  of  200-300  keV/m. 

By  adjusting  the  boundary  values  at  £  in  this  way,  we  can  account  for 

the  slow  variation  of  y_  and  e  with  z  due  to  single  particle  processes, 

L 

while  including  the  fast  variations  due  to  collective  fields  directly 
in  the  beam  equations. 

It  is  frequently  convenient  to  consider  a  "serai-infinite" 


beam,  i.e.  one  which  is  long  enough  so  that  the  plasma  return  current 


has  effectively  decayed  away  at  the  tail,  In  this  case 

UL  =  U»  =  VC' 


(58) 


and  the  boundary  conditions  (53a)  and  (54a)  become  simply 


(52) 


dR 

dC 


dR 

dC 


0. 


(53) 


A  highly  relativistic  beam  of  finite  length  can  always  be  regarded  as 
equivalent  to  the  front  end  of  a  "semi-infinite"  beam  with  properly 
chosen  values,  y^  and  R^ ,  at  C  =  00 ,  since  information  propagates  only 
backward  in  the  beam.  Thus  the  idea  of  a  semi-infinite  beam  is  quite 
general,  and  it  is  convenient  to  use  R  ,  y  and  U  as  scaling  factors 
for  dimensionless  coordinates. 

Although  the  boundary  conditions  (49)  -  (54)  are  sufficient 

to  determine  the  five  quantities  R,  y,  U,  £  ,  and  k  as  functions 

z 

of  £ ,  one  further  constraint  is  needed  to  determine  the  pinch  point 
velocity  8c.  This  additional  condition  follows  from  the  requirement 
that  a  steady  state  exist.  After  the  transient  initial  beam  head 
expansion,  the  beam  axial  profile  should  approach  an  isomorphous  limit 
as  the  radius  at  the  head  RQ  continues  to  increase,  and  8  should 
approach  some  limiting  value.  In  our  time-independent  formulation,  we 
must  therefore  find  solutions  with 


r  ~ 


and  require  that  8  be  insensitive  to  the  exact  value  of  Rq  this  limit. 
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F.  Solution  of  the  Equations 

The  nonlinearity  of  the  beam  equations  and  the  mixed  boundary 
conditions  preclude  solution  either  by  standard  finite-difference 
methods  or  by  numerical  integration.  Instead,  we  use  the 
iterative  method  summarized  in  Fig.  1,  which  is  convenient  for  the 
"semi-infinite"  beam  model  discussed  in  Sec.  HE.  We  define  £  = 

C  -  and  choose  a  test  value  of  $,  which  we  denote  by  B.  We  then 
calculate  the  beam  axial  profile  over  an  interval  0  <  r,  <  L  that  is 
chosen  large  enough  that  all  quantities  effectively  reach  their  £  =  00 
values.  To  start  the  iteration,  we  take  U  to  have  the  form 

U(0)  =  Um  (1  -  exp  (-  5|)  +  |  exp  (-£)]  ,  (59) 

where  we  use  parenthesized  superscripts  to  identify  the  successive 
iterates.  This  expression  satisfies  the  boundary  conditions  UQ  =  0 
and  u  =*  0  ,  and  the  coefficient  5  can  be  chosen  to  give  a  suitable 

L  00 

scale  length.  We  calculate  R  ^  and  by  numerically  integrating 

a  (0) 

the  envelope  equations  (13)  and  (14)  from  ?  =  L,  and  R  is  then  used 

A 

in  integrating  the  conductivity  equation  (43)  from  £  =  0.  since  the  U 
equation  (35)  is  a  linear  ordinary  differential  equation  with  known 
boundary  values,  it  is  readily  reduced  to  tridiagonal  matrix  form  and 
solved  numerically  for  the  first  iterate  .  This  sequence  of  calcu¬ 
lations  is  repeated  until  no  significant  change  in  the  beam  axial 
profile  occurs  with  successive  iterations.  Since  the  only  undetermined 
boundary  value  is  the  beam  head  radius  R^,  a  convenient  criterion  for 
iteration  convergence  is  that  the  quantity  |  R^n+^  -  R^I/R^ 
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be  less  than  some  prescribed  error  limit  for  one  or  more  successive 

iterations.  This  procedure  is  found  to  converge  reliably  for  ic^R^  >  1 

and  R  /R  s  50,  provided  that  the  finite-difference  grid  in  the  U 
0  oo 

calculation  is  fine  enough  to  resolve  the  beam  head  structure. 

To  determine  the  pinch  point  velocity,  we  systematically  vary 

0  to  find  solutions  of  the  beam  equations  with  R^/R^  >>  1.  Figure  2 

shows  a  typical  plot  of  R^  vs  0 .  The  increase  of  R^  with  0  occurs 

because  the  velocity  of  beam  electrons  relative  to  the  pinch  point, 

v  -  V  «  (1-0) c,  is  reduced  as  0  increases.  Consequently,  beam 
z  p 

electrons  traversing  the  pinch  region  are  subjected  to  the  induced  E 

z 

there  for  a  longer  time,  and  the  resulting  reduction  in  y  drives  radial 

expansion.  The  upper  limit  of  0  is  reached  when  the  energy  loss 

2 

approximately  equals  the  initial  kinetic  energy  (y^-Dmc  .  This  0 

value  maximizes  Rq/R^ ,  and  the  corresponding  value  of  0c  should  equal 

the  physical  pinch  point  velocity  =  $c,  provided  that  the  assumptions 

1-0  <<  1  and  R?/r£  «  1  are  also  satisfied.  Because  R„  increases 
c  0  0 

sharply  as  0  0,  the  limiting  value  is  easily  estimated  and,  as 

required,  is  insensitive  to  the  calculated  beam  head  radius.  The 
erosion  rate,  defined  as  the  speed  at  which  the  pinch  point  recedes 


from  the  beam  head,  is  then  v  -  V  ss  (1-6) c. 

z  p 
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III.  Analytic  and  Numerical  Results 

The  basic  equations  (13)  for  R,  (14)  for  Y,  (35)  -  (39)  for  U 

and  $  ,  and  (43)  -  (47)  for  k  are  quite  complicated,  and  the  depend- 
z 

ence  of  the  erosion  rate  and  other  quantities  on  the  beam  parameters  is 

not  immediately  evident  from  these  equations.  However, a  simple  formula 

for  the  erosion  rate  can  be  derived  when  corrections  of  order  X  ^  = 

[£n(R?/R?) I-1  are  neglected.  Using  this  formula  and  again  neglecting 
c 

terms  of  order  X  \  the  equations  for  R,  7,  U,  and  <  can  be  written 
in  a  much  more  concise  dimensionless  form  that  displays  the  parametric 
dependences.  These  simplified  equations,  together  with  data  obtained 
from  numerical  solutions  of  the  equations,  show  that  the  essential 
features  of  the  axial  profiles  R(?)  ,  y(?),  U(c)  and  <(£)  can  be 
described  in  terms  of  scaling  with  the  dimensionless  parameters  which 
control  the  background  conductivity.  This  program  will  be  carried  out 
in  the  remainder  of  this  section.  The  approximate  analytic  results 
and  scaling  relations  will  be  compared  to  extensive  exact  numerical 
solutions  of  (13),  (14),  (35) -.(39),  and  (43) -(47)  obtained  by  methods 


described  in  Sec.  II. 


A.  Erosion  Rate 


When  a  steady  state  exists,  a  scaling  law  for  the  erosion 

rate  follows  almost  immediately  from  the  equations  that  determine 

energy  flow  along  the  beam.  Approximate  scaling  laws  have  previously 
12 

been  obtained  '  by  a  similar  argument ,  in  which  the  ohmic  energy  loss 
of  the  beam  is  equated  with  some  estimated  fraction  of  the  kinetic 
energy  lost  through  erosion. 

The  electron  energy  as  a  function  of  5  is  given  by  the 
integral  of  (14) , 


Y(C)  -  Y, 


0  =  (1-6) me2  j  dC  & 


(5  ) 


(60) 


Here,  the  time  variation  of  y  is  due  to  ohmic  loss  and  the  steady  state 

assumption  has  been  used  to  express  it  as  a  (  variation.  From  (.36)  , 

£  is  given  by 
z 

£„  *  -  fr  (XU)  ,  (61) 

z  d£ 

where  the  last  two  terms  of  (36)  are  neglected,  since  they  are  smaller 

than  the  first  by  a  factor  of  order  X_1<<1.  Substituting  (61)  into 

(60)  and  using  U  ■  0  then  gives  the  dimensionless  erosion  rate 
o 

e  =  e  L«L>M  (62) 

1  B  Y(S)-Y„ 

o 

for  any  S.  If  t,  is  well  behind  the  pinch  point,  then  X(C)  ss  and 

U(£)  as  i  (O/c,  where  the  net  current  I  is  the  beam  current  less  the 
n  n 

plasma  return  current.  Moreover,  y(C)  y>  YQ  for  all  cases  with 
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Y(?)  >>  1/  as  discussed  in  Sec.  IIF  and  shown  in  what  follows.  In 


such  cases ,  (62)  reduces  to 


1-8* 


\Jl  (?)| 

00  n 


(63) 


racJ  y(?)  ' 

which  shows  that  |lJ/Y  is  approximately  constant  behind  the  pinch 
point.  For  a  semi-infinite  beam,  defined  as  one  that  is  long  compared 


with  the  magnetic  decay  length,  we  can  evaluate  (63)  at  z  -  ®,  where 
rn  =  anc*  Y  is  also  known.  We  then  obtain  an  explicit  formula  for 
the  erosion  rate, 


1-8 


x  K 

oo 1  jj' 


(64) 


The  last  expression  relates  the  erosion  rate  to  the  familiar  ratio  of 

I.  to  the  AlfvSn  current  I. . 
b  A 

The  argument  leading  to  (64)  is  valid  only  when  1-8  is  small 

because  the  £  relation  (61)  was  derived  by  assuming  that  0  ss  1.  A 
z 

more  general  derivation  given  in  Appendix  A  shows  that 


1  -  8  « 


(65) 


which  is  valid  for  any  8,  and  reduces  to  (64)  when  1  -  8  <<  1. 

Numerical  solutions  of  (13),  (14),  (35)  -  (39),  and  (43)-(47), 
in  which  the  beam  current  and  ionization  coefficient  A^  are  varied 
over  two  orders  of  magnitude,  substantiate  the  accuracy  of  (64) .  The 
erosion  rate  is  quite  independent  of  the  coefficients  that  govern 
conductivity  <(£),  even  though  these  coefficients  strongly  affect  the 
?  dependence  of  R,  k,  U,  and  y.  The  erosion  rate  is  also  insensitive 
to  the  choice  of  beam  profile  (Gaussian,  Bennett  or  square)  and  to  the 
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parameter  a  in  (33a) ,  which  strongly  varies  the  width  of  the  conducting 
channel.  Table  2  compares  calculated  erosion  rates  with  values  given 
by  Eq.  (64)  for  a  very  broad  range  of  y,  with  =  0  (initially 
neutral  gas)  and  conductivity  coefficients  given  by  Table  1  for  air. 
Other  beam  parameters  are  chosen  to  typify  the  new  generation  of 
linear  induction  accelerators  being  developed  at  Lawrence  Livermore 
Laboratory.  The  calculated  values  of  1  -  8  are  consistently  about 
13%  below  the  estimates.  This  inaccuracy  is  within  the  error  expected 
from  the  estimate  (61)  of  The  scaling  of  1  -  8  with  1^,  y^  and 

X  is  accurately  predicted  by  (64) . 

The  erosion  rate  expression  indicates  that  the  existence 
of  a  steady  non-zero  erosion  rate  depends  on  ohmic  energy  loss:  if 
Y(C)  -  yq  vanishes,  then  (62)  gives  the  physically  meaningless  result 
1  -  8  =  “.  We  have  verified  this  dependence  by  redoing  the  numerical 
solutions  with  y  held  constant,  rather  than  being  determined  by  (14). 
The  steady  state  model  fails  to  produce  a  meaningful  solution  in  this 
case.  We  believe  that  expansion  and  erosion  of  the  beam  head  will 
occur  even  if  y  is  held  constant,  but  not  as  a  steady-state  process. 


28 


B.  Dimensionless  Equations 

Numerical  solutions  of  the  beam  equations  (13),  (14),  (35) — 
(39),  and  (43) -(47)  over  a  wide  range  of  parameters  indicate  several 
approximations  that  can  be  made  with  little  loss  of  accuracy  in  the 
calculated  quantities  R(5) ,  <(?),  U(S),  Y(C)»  and  the  erosion  rate  1-8: 

(i)  The  choice  of  radial  beam  density  profile  n^fr)  has 
little  effect  on  the  results.  When  profiles  are  calculated  using 
Gaussian,  Bennett,  and  square  density  profiles,  the  corresponding 
quantities  typically  differ  by  less  than  3%. 

(ii)  Varying  the  parameter  o  in  (33a) ,  which  controls  the 
shape  of  the  radial  conductivity  profile,  has  a  comparably  small  effect 
on  results,  even  though  a  =  1  corresponds  to  a  conducting  channel  of 
the  same  radius  as  the  beam,  while  a  =  0  corresponds  to  a  constant 
conductivity  out  to  a  radius  Rc  which  is  orders  of  magnitude  larger. 

(iii)  The  terms  in  the  field  equation  (35)  which  arise  from 
displacement  current  terms  in  (15)  and  (16)  also  affect  the  results 
by  only  a  few  percent,  provided  that  A^| X^|  s  0.1.  Neglecting  these 
terms  eliminates  the  second  derivative  in  (35)  and  greatly  simplifies 
the  equation.  We  could  not  have  predicted  the  validity  of  this  approxi¬ 
mation  from  the  equations.  Due  to  the  very  rapid  variation  with  C  in 
the  beam  head,  the  displacement  current  terms  are  not  obviously 
negligible . 

(iv)  Calculated  axial  profiles  are  insensitive  to  parameter 


changes  that  alter  the  erosion  rate  but  leave  the  conductivity 


generation  rate  unchanged,  when  y„  is  varied  between  50  and  1000 
while  other  parameters  are  held  constant,  R(C),  U(C)  and  ie(C)  are 
effectively  unchanged,  even  though  1-6  changes  by  more  than  an  order 
of  magnitude. 

We  now  use  these  results  to  write  the  beam  equations  in  a 
much  simplified  dimensionless  form.  For  specificity,  we  consider  a 
Gaussian  beam  profile  and  a  flat  conductivity  profile  (a  *  0) . 
According  to  (i)  and  (ii) ,  the  results  are  not  sensitive  to  these 
choices.  In  (35)  and  (37)  we  neglect  the  displacement  current  terms 
discussed  in  (iii) ,  as  well  as  terms  that  are  corrections  of  order 
X  1.  We  explicitly  use  (64)  for  1-6  in  Eqs.  (13)  and  (14).  Finally, 
we  choose  a  suitable  set  of  dimensionless  variables. 


0  =  U/U^  ico/^  , 

(66a) 

R  =  V*  , 

OR 

(66b) 

-<  ! 

Ill 

> 

8 

(66c) 

K  =  KR", 

(66d) 

5  SC/V 

(66e) 

Equations  (13),  (14),  (35)-(37),  and  (43)  can  then  be  recast  as 


(67) 


(68) 
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where 


P 


(71) 


In  this  dimensionless  form,  the  quantities  5(0,  y (£),  R (?), 
and  k(C)  depend  only  the  dimensionless  parameters  u ,  ,  1-6, 

and  the  dimensionless  pre-existing  gas  conductivity  =  k^R^  ,  if 
is  nonzero.  However,  the  results  are  insensitive  to  1-6,  as  discussed 
in  (iv) .  This  results  from  the  structure  of  (68) •  Since  R  »  1  and 
5*0  near  the  beam  head,  (68)  reduces  there  to  ydR/dt  »  constant.  In 
the  pinched  region  ydR/dC  is  small  and  slowly  varying,  so  that  the 
envelope  equation  gives  the  quasi-equilibrium  condition  R  «  (yU)1^2. 
For  any  UGT) ,  therefore,  the  erosion  rate  enters  (68)  only  in  the 
smooth  transition  between  these  limits.  Also,  since  we  are  interested 
principally  in  the  beam  structure  around  the  pinch  point,  where 
recombination  is  typically  unimportant,  the  recombination  coefficient 
A 2  can  be  neglected  in  analyzing  beam  scaling.  We  thus  conclude  that 
the  axial  profiles  depend  significantly  on  only  two  parameters,  p 
and  initial  k  value,  k  . 
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C.  Axial  Profile  of  the  Beam  and  Conducting  Channel  for  k 

1.  General  Considerations 

We  have  solved  the  equations  of  Sec.  II  for  cases  ranging 
from  v  =  0.1  to  y  *  100.  The  5  dependence  of  the  principal  beam  and 
channel  features  is  shown  in  Fig.  3-5  for  three  typical  cases  with 
k  =  0  and  y  =  0.32,  3.2,  and  32,  corresponding  approximately  to  1,  10 
and  100  kA  beams  in  initially  un-ionized  air  at  standard  density.  In 
each  case  the  beam  has  two  distinct  regions,  an  expanded  head  and  a 
self-pinched  body,  with  quite  different  structure  and  scale  lengths, 
we  shall  proceed  by  showing  how  each  of  several  beam  properties  differs 
qualitatively  between  the  two  regions  and  defining  a  transition  point 
between  the  two  regions  for  each  property.  These  various  transition 
points  typically  lie  close  together,  as  seen  in  Table  3. 

At  the  beam  head,  where  the  gas  conductivity  is  zero,  the 
magnetic  pinch  force  is  canceled  by  electrostatic  repulsion  due  to  the 
beam  space  charge,  leaving  U  =  0.  As  we  move  back  from  the  head  of  the 
beam,  k  increases  and  the  space  charge  begins  to  be  neutralized  by  the 
radial  outflow  of  plasma  electrons.  Nevertheless,  U  is  so  weak  through¬ 
out  the  head  region  that  the  beam  head  expands  freely,  constrained  only  by 
its  own  inertia.  Within  our  steady  state  model,  this  rapid  beam  head 
expansion  appears  as  a  large  value  of  | dR/d? | .  in  the  self-pinched 
region,  on  the  other  hand,  the  space  charge  has  been  completely 
neutralized,  the  beam  is  in  pressure  balance  with  the  pinch  field  U, 
and  |dR/d?|  is  small.  It  is  therefore  natural  to  define  the  point 
where  |dR/d£|-1  (d2R/d?2)  is  maximal  as  the  pinch  point  ?  . 
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j 

The  rapid  increase  in  U  near  Z  causes  a  substantial 

P 

induced  field,  as  is  evident  from  Eq.  (61) .  Beam  electrons  passing 

through  this  Ez  "spike"  lose  energy  rapidly,  and  conservation  of 

relativistic  momentum  ymr  leads  to  rapid  expansion  of  the  beam  ahead 

of  the  E  spike.  Thus  the  point  ?  at  which  E  is  maximal  typically 
z  e  z 

is  close  to  Z  . 

P 

Since  <  increases  rapidly  with  C ,  the  Ez  spike  drives  a 

plasma  return  current  1^.  Figs.  3-5  show  that  1^ =  I^/I^  increases  rapidly 

just  as  the  space  charge  is  being  neutralized  and  reaches  a  maximum 

at  a  point  which  is  close  to  but  somewhat  larger  than  Z^  because 

of  the  steadily  increasing  conductivity.  As  Z  increases  further  and 

we  enter  the  beam  body,  I  falls  off,  but  if  y  >  1  this  decrease  is 

P 

slow  because  of  the  rapidly  increasing  conductivity.  In  the  beam  body, 

where  there  is  no  space  charge,  the  normalized  pinch  force  U  reduces 

to  c  <rB„>/I,  and  is  equal  to  the  normalized  net  current  I  =  1  -  I  . 

u  d  n  p 

-  1  — o 

U  then  approaches  unity  with  a  scale  length  —  XkR  determined  by 

O 

resistive  diffusion,  as  is  evident  from  (67).  By  contrast,  the  scale 
length  in  the  expanded  head  is  determined  by  the  conductivity  gener¬ 
ation  rate  yR~2  and  is  typically  much  shorter. 

Even  though  Figs.  3c,  4c,  and  5c  show  U(S)  increasing 
smoothly  and  monotonically ,  comparison  of  the  U(5)  and  In(C)  curves 

shows  the  abrupt  change  in  the  nature  of  the  pinch  force  near  Z  :  in 

P 

the  expanded  head,  electrostatic  forces  reduce  U  by  partially 
cancelling  the  magnetic  pinch  force,  but  U  is  purely  magnetic  in  the 
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beam  body.  The  onset  of  plasma  return  current  and  the  rapid  decrease 
of  space  charge  occur  in  the  region  where  k2R2  lies  roughly  between 
8A-1  and  8,  in  which  the  k2R2  terms  in  (67)  begin  to  dominate.  There¬ 
fore  the  point  where  Atc2R2/8  =  1  represents  another  transition  point 

for  the  electromagnetic  fields,  and  is  typically  also  close  to  C  . 

P 

Since  the  pinch  force  is  in  quasistatic  equilibrium  with 

the  beam  pressure  throughout  the  beam  body,  (68)  gives  UyR2  Ss  1  in 

this  region.  In  the  expanded  head,  where  the  pinch  force  is  too  weak 

to  cpnfine  the  beam,  UyR2  <<1.  Figures  3b,  4b  and  5b  show  a  fairly 

sharp  transition  between  these  regimes.  We  arbitrarily  define  a 

transition  point  ?  where  UyR2  =  0.95.  The  location  of  ?  is  of 
qs  qs 

some  significance,  since  it  delineates  the  regime  in  which  a  quasi¬ 
static  model  of  beam  equilibrium  is  valid,  for  use  in  stability  theories, 

for  example.  Although  5  is  always  fairly  close  to  C  ,  we  find  that 

qs  p 

C  „  <  S  for  y  s  1  and  C  >  £  for  y  >  1. 
qs  p  qs  p 

The  relative  value  of  (  ,  j  (  and  c  are  given  in 

p  n  i  e  qs 

Table  3 ,  over  a  large  range  of  y .  Several  other  beam  properties  are 
also  summarized  in  Table  3,  and  will  be  discussed  in  the  next  sections. 
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2 .  Validity  of  the  Steady  State  Assumption 

In  an  actual  beam,  the  radius  R  at  the  beam  head 

o 

increases  constantly  with  time,  which  is  not  strictly  consistent  with 

our  steady  state  treatment.  However,  our  treatment  is  based  on  the 

premise  that,  when  Rq  has  become  much  larger  than  the  pinched  radius, 

the  value  of  Rq  has  a  negligibly  small  effect  on  all  beam  properties 

in  the  region  of  interest,  where  R  <<  R  .  This  premise  was  supported 

o 

by  the  observation  in  Sec.  IIP  that  as  8  ->■  B,  large  changes  in  Rq 

correspond  to  very  small  changes  in  6.  It  is  also  necessary  that  R(5), 

y(C)  and  <(£)  converge  to  steady  state  values  when  R  >  R  ,  where  R 

op  p 

is  the  radius  at  the  pinch  point.  This  requirement  is  also  found  to  be 

well  satisfied  in  the  numerical  solutions.  As  an  example,  Fig.  6  shows 

R(S/RJ  for  a  typical  case,  obtained  by  the  method  described  in  Sec. IIP, 

using  three  slightly  differing  value  of  S .  Even  though  the  8  variation 

results  in  very  different  values  of  R  ,  R(5/R  )  profiles  are  nearly 

o  00 

identical  when  R^  >  10R^.  We  may  therefore  describe  the  time  during 
which  the  beam  propagates  from  its  source  until  Rq  reaches  the  order 
of  10Rp  as  a  period  of  transient  beam  head  expansion.  Time-dependent 
treatments  have  indicated  that  this  initial  transient  expansion  occurs 
considerably  faster  than  the  steady  state  erosion  that  follows.^- 
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In  the  beam  body,  the  degree  of  beam  pinching,  i.e.  the 

value  of  the  dimensionless  radius  R,  depends  on  the  fractional  plasma 

return  current  I  .  Although  it  is  always  true  that  R  -*•  1  and  i  -*•  0 
P  P 

as  5  becomes  very  large,  the  values  of  R  at  ?  and  of  I  at  its 

p  p 

maximum  point  are  indicative  of  conditions  over  a  substantial 
length  of  the  beam  body:  A  beam  is  tightly  pinched  throughout  the  body 
if  R(?  )ssl,  whereas  a  beam  with  R(S^)  >>  1  may  be  described  as  weakly 
pinched.  The  numerical  results  tabulated  in  Table  3  indicate  that 

R(?  )  «  1  +  y  (72) 

P 

for  0.1  v<  M  100.  A  beam  is  therefore  seen  to  be  well-pinched  if 
y  <  1  and  weakly  pinched  if  y  >  1,  in  agreement  with  many  experimental 
observations  that  beams  tend  to  be  weakly  pinched  when  the  conductivity 
rises  too  rapidly. 

Figures  3-5  and  Table  3  further  suggest  that  the  axial 

profiles  R(S),  U(?)  and  k  (? )  can  be  characterized  in  the  region 

£  S  ?  by  very  approximate  self-similar  expressions  with  y  as  a 
P 

scaling  parameter.  To  obtain  such  scalings,  we  assume  forms 
n 

R  =  y  f1(ymC)  ,  (73a) 

n 

5  =  y  f2(ym?)  ,  (73b) 

n 

<  -  y  f3(ym?)  ,  (73c) 
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and  we  neglect  both  the  recombination  term  A2p2k2  in  (70)  and  the 

k2 R2  terms  in  (67).  We  choose  the  coefficient  n^  to  be  consistent  with 

(72)  and  require  that  n2  =  ”  n  ?  so  that  the  beam  can  reach  quasistatic 

equilibrium  at  ?  «  To  satisfy  the  requirements  of  self -similarity, 

we  must  also  choose  the  radius  R  at  the  beam  head  to  be  some  fixed 

o 

multiple  of  R^.  The  only  possible  scaling  of  (67)  -  (70)  that  is  con¬ 
sistent  with  these  requirements  is 


R=f1(y1/2C),  (74a) 

5  =  (74b) 

k  *  y1/2f3(y1^20  ,  (74c) 

for  u  «  1,  and 

R*pF1(?),  (75a) 

U  =  u_1F2(0,  (75b) 

<  =  U_1F3  (?)  ,  (75c) 


for  y  »  1.  Neither  scaling  is  exact,  but  both  are  useful,  and  the 
numerical  solutions  illustrated  in  Figs.  3-5  and  Table  3  verify  their 
approximate  validity. 

These  scaling  laws  yield  several  types  of  information  about 
the  beam  evolution,  (i)  They  indicate  how  the  axial  profile  of  the 
beam  depends  on  the  parameter  y :  the  "trumpet"  shape  of  the  expanded 
head  becomes  more  extreme  as  y  increases.  (ii)  We  have  argued  that  the 
beam  passes  through  a  transient  period  of  very  rapid  expansion,  during 
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which  R  increases  to  approximately  10R(?  ) ,  before  the  steady  state  is 
o  p 

established.  The  extent  in  z  of  the  expanded  region  is  thus  equal  to 

C  when  R  'v  10R(5  ).  The  scaling  relations  (74)  and  (75)  indicate 

pop 

that  this  region  scales  as  y1/2  for  y  <  1  and  reaches  a  constant  value 

for  y  >  1.  However,  (75)  also  indicates  that  the  body  of  the  beam 

expands  a  great  deal  if  y  »  1.  (iii)  An  electron  which  was  originally 

located  in  the  pinched  body  of  the  beam,  ?  >  ?  ,  but  which  later  finds 

P 

itself  in  the  expanding  nose  as  the  pinch  point  moves  back  into  the 
beam,  will  in  general  escape  radially  at  a  velocity  less  than  its 
original  transverse  thermal  velocity,  because  the  expanding  beam  does 
work  against  the  magnetic  field.  According  to  (74)  and  (75) ,  the 
escape  velocity  scales  as 

g-  a-«le|§l  u  <  1  (76a) 

«  yXjlJ/X,,  y  >  1.  (76b) 
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.  Axial  Profiles  and  Scaling  for  *  0 

When  the  background  gas  is  sufficiently  preionized  to 

neutralize  the  beam  space  charge,  the  bean  still  has  a  freely  expanding 

head  and  approaches  steady-state  erosion.  However,  the  fields  are 

quite  different  from  a  case  with  *  0,  and  we  find  significantly 

different  scaling  laws.  Since  Ak2R2/8  s  1  for  £  <  £  when  c  -  0,  a 

p  o 

condition  that  k  dominate  beam-generated  conductivity  for  all  C  <  C 
o  p 

is  that 

\i<  2R2  /8  >  1 ,  (77a) 

o  po 

where  R  is  R(C  )  for  the  beam  when  k  =0.  Using  (72)  for  R  ,  this 
po  p  o  po 

condition  becomes 


Ak  2/8  >  1  ,  if  v  <  1,  (77b) 

o 

Ak  2y2/8  >  1  ,  if  u  >  1.  (77c) 

o 

When  (77)  is  satisfied,  the  beam  is  space-charge  neutralized 

essentially  everywhere ,  so  that  U  =  I  .  Nevertheless,  U  =  0  at  the 

n 

beam  head  because  the  plasma  return  current  generated  there  completely 
cancels  the  beam  current.  The  absence  of  a  pinch  force  at  the  beam 
head  leads  to  free  expansion  of  the  head  and  subsequently  to  steady 
state  erosion,  just  as  in  the  case  with  =  0.  Although  the  erosion 
rate  is  independent  of  both  kq  and  y,  as  we  have  seen  in  Sec.  IIIB, 
the  R(5)  profile  is  modified  by  icq.  As  long  as  k  w  kq,  the  scale 
length  for  the  rise  of  In  is  ^  Ak^R2.  The  only  possible  self-similar 
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(78a) 


scaling  for  R(?)  ,  5  (?),  and  U(?)  is  then 
n 

R(C)  =  f.  (?/k  ) 

1  o 

I  (?)  =  U(?)  =  f,(?A  ),  (78b) 

n  2o 

where  f  and  f  2  are  independent  of  y .  The  beam  pinches  down  more 

slowly  than  in  un-ionized  gas,  and  so  for  any  value  of  ?,  R(?)  is 

larger  than  it  would  be  for  k  *  0.  These  features  are  illustrated  by 

o 

the  plots  of  R(?)  and  J  (?)  in  Fig.  7  for  cases  with  y  =  0.2  and  in 

n 

Fig.  8  for  y  =  5. 

For  cases  in  which  is  nonzero  but  not  large  enough  to 

satisfy  (77) ,  the  pre-existing  conductivity  can  modify  the  beam 

properties  in  a  limited  region  of  the  expanded  head,  but  will  have 

little  or  no  influence  on  the  beam  body.  For  sufficiently  large  R  , 

o 

even  a  small  can  result  in  a  value  of  kr  near  the  beam  head  large 

enough  to  neutralize  the  space  charge  there.  However, 

decreases  at  large  £  values  as  R  decreases,  and  space  charge  forces 

can  reappear  at  some  ?  <  ?^.  This  results  in  the  somewhat  complicated 

I  (?)  curves  for  <  *  0.07  in  Figs.  7b  and  8b.  When  y  <  1,  Eq.  (74) 

n  o 

indicates  that  the  beam  tends  to  pinch  down  faster  when  conductivity 

is  increased.  Thus  a  small  non-zero  value  of  increases  | dR/d? | 

near  the  beam  head,  just  as  a  small  increase  in  y  would.  However,  in 

the  part  of  the  beam  where  ic  >  ic  the  R(?)  curves  have  the  form  they 

o 

would  have  for  <  *  0  but  are  shifted  toward  the  head,  as  Fig.  7a  shows, 

o 

In  particular,  R(?)  at  and  behind  the  pinch  point  is  unchanged  by  <  . 
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When  |i  >  1,  dR/dC  is  insensitive  to  small  changes  in  conductivity, 
and  if  (77)  is  not  satisfied  the  entire  R(C)  curve  is  essentially 
independent  of  as  seen  in  Fig.  8a. 
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IV.  Summary 

The  pinch  £rame  formalism  presented  here  is  an  efficient  and 
conceptually  simple  method  for  calculating  the  beam  head  erosion  rates 
and  axial  profiles  of  relativistic  electron  beams.  The  principal 
shortcomings  of  the  method  are  the  non-selfconsistent  treatment  of 
emittance  and  energy  changes  due  to  collisions  and  radiation*  and  an 
inability  to  model  early  transient  expansion  of  the  beam  head.  None¬ 
theless  the  formalism  gives  usable  information  about  beam  structure 
and  an  unambiguous  erosion  rate  value  without  requiring  time-dependent 
calculations . 

We  have  examined  relativistic  electron  beams  with  a  wide  range  of 
parameters  and  can  make  several  generalizations: 

1)  If  energy  loss  and  expansion  from  collisions  and  radiation  are 
neglected,  erosion  continues  at  a  steady  rate  after  initial  transient 
expansion  of  the  head.  For  a  constant-current  pulse,  this  rate  is 

lyol 

y7F5 - 17k &  c, where  5  is  any  point  well  behind  the  pinch  point.  The 

radial  dependences  of  beam  density  and  conductivity  and  the  maximum 
radius  of  the  conductivity  channel  Rc  affect  erosion  rates  only  weakly 
by  altering  the  beam  inductance  A  . 

2)  when  the  background  gas  is  initially  un- ionized  and  the 
current  is  high  enough  that  tc^R^  >>  1,  then  the  beam  axial  profile 
near  the  beam  head  is  controlled  principally  by  the  dimensionless 
parameter  u  =  47rR2c  ^d<a>/d?  that  controls  conductivity  generation. 

For  air  with  ionization  caused  by  beam-molecule  collisions, 
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U  as  3  x  10  IjJ  (2eV/Te>1/2.  When  y  <<  1,  the  peak  return  current  is 

small  compared  with  1^,  and  the  beam  pinches  to  R  as  R^.  In  contrast , 

the  return  current  is  comparable  with  1^  when  y  >>  1  and  the  radius 

at  the  pinch  point  is  approximately  yR^. 

3)  When  the  initial  gas  conductivity  k  is  large  compared  with 

o 

beam-generated  conductivity  throughout  the  beam  head,  the  beam  space 

charge  and  current  cure  neutralized  essentially  at  The  beam  pinches 

with  a  scale  length  Xk  B?/Q.  Although  there  is  no  distinct  pinch 

o 

point,  beam  head  erosion  proceeds  at  the  same  rate  as  for  <  =  0. 
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appendix 


The  erosion  rate  calculations  in  Section  II IA  are  limited  to  con¬ 
ditions  for  which  (JL-g)  <  0.2  by  the  Bssl  approximation  made  in  deriving 
the  field  equations  (35) -(39).  We  can  estimate  the  rate  for  lower 
values  bY  using  the  field  equations  (23) -(25),  which  do  not 
assume  Bssl.  Integrating  (24)  over  the  entire  beam  gives 
00 

r  dc  E2(r,?)  =|  [Bljr)  +  (1-B)  IQ(r)]  ,  («l) 

5 

where  we  have  used  the  boundary  values  E  £r,0)  =  E  (r,“)  =0  and 

z  z 

u(r,«)  =  rBQ (r,»)  and  have  eliminated  rBg  with  (25).  If  we  assume  that 
the  left  side  of  (Al)  has  the  same  radial  dependence  as  I(r),  then  the 
equation  can  be  radially  averaged  at  ?  *  «  to  give 

00 

■  ajr’T  J «'  ipfrT r  -  |  *  u-B><io>]  <«> 

C 

where  <(  1^  ( r )  /  =  y  1^  for  any  of  the  radial  density  profiles  treated 

here.  Since  IQ(r)  is  the  current  within  radius  r  at  the  expanded  beam 

head,  the  average  of  I  over  the  number  density  at  5  ■  00  is  smaller 

than  1/2  by  approximately  the  ratio  R2/R2  <<  1,  and  the  <1  >  term  in 
D  oo  O  O 

(A2)  is  negligible.  If  we  again  ignore  the  change  of  inductance  along 
the  length  of  the  beam,  then  we  may  set  A(C)  «  X  in  (A2)  and  obtain 

OO 
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;  bx  ixj 

J  «'  (e;>  *  — 


(A3) 


Substituting  this  approximation  into  (3.1)  gives 

•U'bl 


Ay 


*  JL 


1-6  mc3 


(A4) 


which  differs  from  the  8  =  1  result  (3.4)  only  by  the  B  factor  on  the 
right  side.  The  corresponding  expression  for  the  limiting  erosion 
rate  is 


This  expression  agrees  with  (64)  for  1-8  <<  1  and  should  be  valid 
over  a  wider  range  of  1 1^1  and  The  principal  restriction  on  (A5) 

is  that  y,,,  must  be  large  enough  that  the  ultrarelativistic  approxi¬ 
mation  v  %  c  is  satisfied  over  most  of  the  beam, 
z 


45 


Acknowledgment  s 

We  are  grateful  to  Adam  T.  Drobot  for  useful  discussions  of 
numerical  methods. 

This  work  was  sponsored  by  the  Defense  Advanced  Research  Projects 
Agency  (DOD)  under  ARPA  Order  No.  3718,  Amendment  No.  12. 


46 


References 


1.  R.  Biegalski  (private  communication) . 

2.  E.  P.  Lee  (private  communication) . 

3.  K.  A.  Brueckner  (private  communication) . 

4.  R.  F.  Hubbard,  D.  S.  Spicer,  J.  C.  Guillory,  and  D.  A.  Tidman 
"Classical  Model  of  High-Current  Beam  Propagation" ,  JAYCOR  Report 
350-79-012  (1979) . 

5.  R.  R.  Johnston,  R.  L.  Feinstein,  D.  E.  Maxwell,  E.  R.  Parkinson, 
and  E.  E.  Simpson  (private  communication) . 

6.  M.  I.  Haftel,  M.  Lampe,  and  J.  B.  Aviles,  Phys.  Fluids  22^  2216 
(1979) . 

7.  E.  P.  Lee  and  R.  K.  Cooper,  Particle  Accelerators  7_,  83  (1976)  . 

8.  W.  H.  Bennett,  Phys.  Rev.  45^  890  (1934). 

9.  B.  Rossi,  High  Energy  Particles  (Prentice-Hall,  New  York,  1952) , 
pp.  27  and  55. 

10.  B.  Rossi,  Qp.Cit. ,  p.67. 


47 


Table  1.  Typical  v  /p  values  and  direct-ionization  coefficients 
m 

for  air.  He,  Ne,  and  Ar.  A  background  electron 
temperature  of  2eV  is  assumed. 

Gas 


Air 

5 

2.0 

X 

io12 

4.3 

X 

io20 

0.24 

12 

20 

He 

5 

1.3 

X 

10 

2.9 

X 

10 

0.25 

Ne 

5 

3.2 

X 

io11 

1.3 

X 

1020 

0.45 

Ar 

5 

2.0 

X 

1011 

2.4 

X 

1019 

0.14 

Beam 

Energy 

(MeV) 


v  /p 

m  -1 
(see  L) 


,  -1  -1, 

(cm  sec  ) 


“1 

2,  -1. 


2.0  x  10 


6.0  x  10 
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Table  2.  Comparison  of  calculated  erosion  rates  with  analytic 
estimates  for  R  =  0.2  cm  and  R.  =  20  cm.  The  para- 
meters  Ib  and  R^  where  chosen  to  fall  in  the  range 
spanned  by  the  new  generation  of  linear  induction 
accelerators  at  Lawrence  Livermore  Laboratory,  the 
Experimental  Test  Accelerator  (y  'v  10)  and  Advanced 
Test  Accelerator  (y  'v  100)  and  to  cover  a  broad  energy 
range. 


Density 

Profile 

Ifa  (kA) 

Calculated 

i-e 

Analytic 

1-B 

Gaussian 

10 

1000 

4.7  x 

10  " 

5.3  x 

io'3 

10 

500 

9.4  x 

IQ"3 

1.1  X 

<N 

1 

O 

i—A 

10 

200 

2.4  x 

io-2 

2.7  x 

io*2 

10 

100 

4.8  x 

10~2 

5.5  x 

io"3 

10 

50 

9.6  x 

io-2 

1.1  X 

io-1 

10 

20 

2.6  x 

io-1 

2.8  x 

io-1 

20 

100 

9.6  x 

O 

1 

ro 

1-1  X 

io-1 

5 

100 

2.4  x 

io-2 

2.7  x 

io-2 

Bennett 

10 

100 

4.3  x 

io-2 

5.5  x 

io-2 

Square 

10 

100 

5.4  x 

io-2 

5.5  x 

io-2 

49 


Table  3.  Profile  data  for  beams  in  initially  un-ionized  gas.  In 
each  case  8  0.05.  For  p  =  0.1  to  32.  Rc  «  10R(Sp)  is 

chosen  to  facilitate  comparison  with  the  self-similar 
scaling  laws  (74)  and  (75).  For  y  =  100,  Rq  <<  10R(Cp) 
is  used  for  reasons  of  numerical  efficiency,  and  as  a 
result  all  ?  points  are  shifted  toward  zero. 


(a)  R  (b)  j 

n 

Fig.  8  —  Variation  of  beam  axial  Profile  with  <  when  y  =  5 
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